High-velocity collision of two black holes 
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We study nonaxisymmetric collision of two black holes (BHs) with a high velocity v = \dx^ /dx^''\ — 
0.6-0. 9c at infinity, where x'^ denotes four-dimensional coordinates. We prepare two boosted BHs for 
the initial condition which is different from that computed by a simple moving-puncture approach. 
By extrapolation of the numerical results, we find that the impact parameter has to be smaller that 
~ 2.5G'Mo/c^ for formation of a BH in the collision for w — > c, where Mo(? is the initial total ADM 
mass energy of the system. For the critical value of the impact parameter, 20-30% of mass energy 
and 60-70% of angular momentum are dissipated by gravitational radiation for v = 0.6-0. 9c. 
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PACS numbers: 04.25.D-, 04.30.-w, 04.40.Dg 

I Introduction: Clarifying formation process of mini 
black hole (BH) in higher-dimensional spacetimes has be- 
come an important issue since a possibility of BH forma- 
tion in accelerators was pointed out. If our space is the 
3-brane in large [l| or warped [2^ extra dimensions, the 
Planck energy could be of O(TeV) that may be accessi- 
ble with planned particle accelerators. In the presence 
of the extra dimensions, a BH of very small mass may 
be produced in the accelerators and the evidence may be 
detected. 

The possible phenomenology of a BH produced in ac- 
celerators was first discussed in jSi] (see [^] for reviews). 
During the high-energy particle collision of sufficiently 
small impact parameter in a higher-dimensional space- 
time, two particles will merge to form a distorted BH, 
and then, it settles down to a quasistationary state af- 
ter emission of gravitational waves. The quasistationary 
BH will be soon evaporated by the Hawking radiation, 
implying that the quantum gravity effects will be impor- 
tant. The evaporation and quantum gravity effects Q 
have been studied for yielding a plausible scenario (cf. 
P for related issues). By contrast, the analyses for BH 
formation and subsequent evolution by gravitational ra- 
diation have not yet been done in detail (but see [3]). 
These phases are described well in the context of general 
relativity Q, but due to its highly nonlinear nature, any 
approximation breaks down. Thus, numerical relativity 
simulation is the unique approach to study this phase. 

In this paper, we present a new numerical-relativity 
study for high-velocity collision of two BHs in four di- 
mensions. This is the first step for understanding the 
high-velocity collision of two BHs in higher-dimensional 
spacetimes. We perform simulations for two equal-mass 
BHs of no spin. A new approach is adopted for preparing 
the initial condition (see Sec. II). The velocity of each 
BH is chosen in the range, 0.6-0. 9c (where c is the speed 
of light), and the results are extrapolated to infer a result 
for V ^ c. We determine the largest value of the impact 
parameter for formation of new BH and approximately 
determine the final mass and spin of the BH. In the fol- 
lowing, we use the geometrical units in which c = G = 1. 

II Initial condition: There are several methods for 
preparing initial condition for high-velocity collision of 



two BHs. A popular method is the moving-puncture ap- 
proach which has been adopted in a recent work for the 
head-on collision of two high- velocity BHs [7|] . In the sim- 
ple moving-puncture approach in which the three-spatial 
hypersurface is assumed to be initially conformally flat, 
the BHs are not in a stationary state in their own comov- 
ing frame and hence a large amount of spurious gravita- 
tional waves are included, as pointed out in [3]. To avoid 
this unsuitable property, in this paper, we use a differ- 
ent approach from the simple moving-puncture one: We 
superimpose two boosted BHs, as described below. 

The line element of a nonrotating BH in the isotropic 
coordinates is written as 



ds^ — —a^dtl -I- V'o(^^o + ^?/o + dz^), 



where 
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rg — Xq -\- Uq -\- Zq and mo is the BH mass in the rest 
frame of the BH. By boosting the BH in the x-axis di- 
rection with the speed v, the line element becomes 

ds2 ^ _r2(ao - il)Qv'^)dt^ + 2T'^v{al - ^^)dtdx 

+ i;^{Bldx^ +dy^ +dz^), (3) 

where the new coordinates x^^ are related to the original 
one by the Lorentz transformation as t = T{to -\- vxq), 
X — T{xq + vto), y — yo, and z = zq. T is the Lorentz 
factor, r = — v^, and 



(4) 



Note that at i = 0, ro = y/T'^x'^ y'^ -\- z'^. 

Because the exact solution of a boosted BH is described 
by Eq. ([3]), in the following, we consider initial data for 
which the spatial metric has the form 



^/{B^dx^ +dy^ + dz'^). 



(5) 



Before going ahead, we summarize the lapse function (a) , 
nonzero-component of the shift vector (/?*), and non-zero 
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components of extrinsic curvature (Kij) of the boosted 
BH at t = 0: 
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The dash (') denotes the ordinary derivative with respect 
to rp (e.g., = dao/dro), and Kij is derived from 
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Di denotes the covariant derivative with respect to three 
metric, 7^, and we use the relation dtjij = —TvdxoJij- 

Now, we describe initial data for two BHs. Although 
we adopt initial condition which approximately satisfies 
constraint equations in general relativity in this paper, a 
general framework is first summarized. 

We write the conformal factor as 



4' = V'main 
V'main = 1 ' 



2^1 



m2 
2^' 



(12) 
(13) 



where denotes mass parameter of each BH, = 
\/r2(x - Xq)2 + (y - yaY + z2, and {xa,ya,Q) denotes 
the location of each BH at t = 0. Namely, we express 
each BH by a moving-puncture framework in a modified 
form, in which 7^ ^/ {x — XaY + (j^ ~ VaY + ■ is a 
correction term which should be determined by solving 
the Hamiltonian constraint. 

This paper focuses on the equal-mass case in which 
mi = rn2 = TOq, X2 = -Xi = Xp (> 0), y2 = -yi = 
h/2 (> 0). Two BHs are assumed to have the same ab- 
solute velocity but move in the opposite directions each 
other; i.e., vi = —V2 = v {> 0). Here, b denotes the im- 
pact parameter. The total mass energy Mq and angular 
momentum J of the system are 

Mo = 2mor and J = moTvb, (14) 

and the nondimensional spin parameter of the system is 

J bv 

(15) 
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It is natural to expect that a new BH is formed after the 
collision whenever J/Mq < 1, i.e., b < AmoT/v. 

Taking into account that the line element of a boosted 
BH is written by Eq. we write as 



B^ = 



2ri 2r2 ' 
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FIG. 1: (a) Violation of the Hamiltonian constraint along x 
axis for xo = 160mo, 6 = 0, and v = 0.9. The violation is 
defined by |//|/(|//i | I//2I -h I//3 1) where Hi = H2 = 

-i?i/)/8, H3 = [KijK'' - (A'fc*)2],/,V8, and H = Hi + H2 + 
Hi should be zero if the Hamiltonian constraint is satisfied. 
Di and R are the covariant derivative and Ricci scalar with 
respect to the conformal three metric ^ij = 'tl)~*'^ij. (b) The 
same as (a) but for numerical results near one of BHs just 
before coUision for xo = 160mo, bv/(moV) = 5.088, and v = 
0.8 in the simulations of different grid resolution, h/mn = 
0.0625 (dotted curve; a,t t = 256mo), 0.050 (dashed curve; 
at t = 256mo), and 0.042 (sofid curve; at f = 258mo). BH 
is located for 3 ^ x/mo < 4. Note that convergence is not 
seen because the degree of the constraint violation seems to 
be primarily determined by the initial condition. 



For the extrinsic curvature, we basically superimpose two 
parts as 



K,, = K 



lij 



2i] 



SK, 



(17) 



where Kaij (a = 1,2) are defined from Eqs. ([7|)- ([T(I)l by 
replacing tq to (a = 1,2). SKij is a correction term 
which should be determined by solving the momentum 
constraint. 

In this paper, we adopt an approximate initial condi- 
tion of (/) = and SKij — 0. The adopted initial data 
does not satisfy the constraint equations of general rela- 
tivity for finite values of Xa and ya or for nonzero value of 
V. However, for the case that Ra = {x-^ + y^y^'^ /mo '> 1, 
the violation of the constraints is tiny because the mag- 
nitude of the violation is proportional to mo/Ra- In the 
present work, we choose Ra/mo ^ 100 (typically 160) 
for the initial condition. In such case, the violation of 
the constraints is ~ mo/Ra = 0(0.01) for most of re- 
gion (see Fig. [1]), and thus, the initial condition approx- 
imately satisfies the constraints. The exception occurs 
around the puncture for which the violation is large, but 
the worst region is hidden inside apparent horizon and 
hence does not play a bad role. The constraint viola- 
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FIG. 2: Area of apparent horizon as a function of time before 
collision for v = 0.9 and bv/(mor) = 4.708 with different 
grid resolution, h/mo ~ 0.0625 (dotted curve), 0.050 (dashed 
curve), and 0.042 (solid curve). 



tion does not disappear during the evolution, but the 
magnitude of the violation remains roughly in the same 
small level at least before collision of two BHs (see Fig. 
[Ijb)). In this method, the BHs are approximately in a 
stationary state in their own comoving frame and hence 
a large amount of spurious gravitational waves are not 
included, in contrast to the initial data prepared by the 
simple moving-puncture approach. 

With this initial data, apparent horizons are located 
approximately for Va = at t = 0. The area of each 

horizon is « IGttTOo within the error of O(10~^) for Xa ^ 
lOOmp (thus the bare mass of the BHs would be « tuq, 
and hence, a fraction of kinetic energy in the total mass 
of the system is 1 — T~^). Furthermore, the area of the 
apparent horizon remains approximately Wttitiq during 
evolution before collision (cf. Fig. [2|). Therefore, it is 
reasonable to expect that numerical solution obtained by 
the simulation provides an approximate solution within 
an error of 1%. In the future work, we will perform 
simulations using improved initial data which satisfies the 
constraints, obtained by computing correction terms, (f> 
and 5Kij, to strictly validate the present strategy. 
Ill Numerical results: For numerical simulation, 
we use SACRA code recently developed by our group 
0. In SACRA, the Einstein equations are solved in a 
modified version of BSSN (Baumgarte-Shapiro-Shibata- 
Nakamura) formalism [10] with a fourth-order finite dif- 
ferencing scheme in space and time and with an adaptive 
mesh refinement algorithm (at refinement boundaries, 
second-order interpolation scheme is partly adopted and 
hence the convergence may reduce to be second order). 
The moving-puncture approach is adopted for following 
moving BHs Gravitational waves are computed 

by extracting the outgoing part of the Newman-Pcnrose 
quantity (the so-called ^'4). Properties of the BHs such 
as mass and spin are determined by analyzing area and 
circumferential radii of apparent horizons. The details of 
our schemes, formulation, gauge conditions, and meth- 
ods for the analysis are described in This reference 
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FIG. 3: Trajectories of relative position of moving punctures 
for V — 0.9 and fe/mo = 6.0, 6.2, and 6.4 (solid, dashed, and 
dotted curves), {bv/mor = 4.708, 4.865, and 5.021.) 



also shows that SACRA can successfully simulate merger 
of two equal-mass BHs. Because an adaptive mesh re- 
finement algorithm is implemented, the moving BHs can 
be computed accurately by preparing a high-resolution 
domain appropriately around the BHs; in the present 
work, we prepare 10 refinement levels. Indeed, we per- 
formed test simulations in which a single high-velocity 
BH is boosted with v = 0.8 and 0.9, and found that 
our code can follow such high- velocity BH for more than 
SOOtoq! e.g., we checked that the area of apparent horizon 
converges at approximately fourth order with improving 
grid resolution and the error is within ~ 0.1% and 1% 
for V — 0.8 and 0.9, respectively, for the grid resolution 
with h = mo/20 (e.g.. Fig. 

Numerical simulations are performed for v — 0.6, 
0.7, 0.8, and 0.9 and xo/mo = 160 changing the im- 
pact parameter b. As expected from Eq. two BHs 
should merge after collision for b < AmoT/v, and hence, 
the critical value of the impact parameter for the BH 
formation (hereafter &crit) should be close to AnioT/v. 
For this reason, the value of b is chosen in the range, 
4 < 6t;/(mor) < 5.5, with the basic step size of b being 
O.lmo (near b = 6crit, the step size of 0.05mo is partly 
used). 

Numerical results depend only weakly on the initial 
separation for given values of v and 6, and for a given grid 
resolution. Indeed, we performed detailed test simula- 
tions for V = 0.8 and 0.9, and xo/mo = 80, 128, and 160. 
If Xo/mo is changed from 160 to 128, the value of 5crit, 
which is one of the most important outputs in this work, 
decreases only by < 0.05mo, and from xo/mo = 160 to 
80, by < O.lmo. Recall that the larger value of the initial 
separation results in the smaller initial constraint viola- 
tion. Thus, the weak dependence of the numerical results 
on xo indicates that the initial constraint violation only 
weakly affects the numerical results. 

The numerical simulations are performed for different 
grid resolutions as h/mo — 0.075, 0.0625, 0.050, and 
0.042. The outer boundaries along each axis are located 
at L « 770mo for all the grid resolution. The value of 
L/c is longer than the duration that the simulation is 
carried out, and hence, spurious effects from the outer 
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FIG. 4: Summary of the final outcomes after the coUision of 
two BHs in the parameter space of {v,b). The circles denote 
that a BH is formed after collision, whereas the crosses do not. 
The solid curve denotes b/(moT/v) — 4.0, and the dashed 
curves are 4.9, 5.0, and 5.1. 



boundaries are excluded. We find that the value of fecrit 
depends only weakly on the grid resolution (it increases 
only by ~ O.lmo 1.5%) if we change h/rriQ from 0.0625 
to 0.042). The area of a BH formed after the merger, 
and total energy and angular momentum dissipated by 
gravitational waves (Ai? and A J) depend more strongly 
on the grid resolution. However, our results indicate a 
convergence slightly better than second order with im- 
proving the grid resolution as long as ft,/mo < 0.0625 
and 0.6 < v < 0.9, although the convergence becomes 
slow for V ^ 0.9. This is natural because the coordinate 
radius of the apparent horizon in the x-axis direction 
for each BH is proportional to T^^. Obviously, the bet- 
ter grid resolution is necessary for the ultra-high velocity 
collision w — > 1, but this is beyond scope of this paper. 

Figure [3] plots trajectories of moving punctures for 
V = 0.9, xo = 160mo, and bv/{mQT) = 4.708, 4.865, 
and 5.021. Here, we plot relative position [(3:2 — 
xi)/2mo, {1/2 — yi)/2mo] on the orbital plane. For the 
first two impact parameters, the BH is formed after the 
collision, whereas for the other, each BH escapes from 
the center after scattering. For the small values of b, two 
BHs form a bound orbit when the orbital separation be- 
comes small enough. For a sufficiently small value of b, 
such as bv/{moT) ^ 4.7, two BHs merge within one orbit. 
For a value of b close to 6crit, two BHs rotate around each 
other for more than one orbits before two BHs merge, as 
shown in Fig. [31 For b > 6ciit, two BHs rotate around 
each other for small separation. However, they do not 
constitute a bound orbit because of the large centrifugal 
force, and eventually, each BH escapes from the center. 

To summarize the outcomes after the collision, we gen- 
erate Fig. m which shows a parameter space composed of 
{v,b). We plot the circles for the case that two BHs 
merge after the collision, whereas the crosses are plot- 
ted when two BHs do not merge. The solid curve de- 
notes b = AmoT/v for which the nondimensional spin 
parameter of the system is unity at i = 0. The three 
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FIG. 5: Outgoing part of the Newman-Penrose quantity for 
V — 0.8 and bv/{moT) — 5.088. The solid and dotted curves 
denote I = m — 2 and I — m = 4 modes for the best-resolved 
run, and the dashed curve denotes I = m = 2 modes for the 
second-best run. tret = approximately corresponds to the 
time of the onset of merger. D denotes a distance between 
the source and observer. 



dashed curves denote bv/moT — 4.9, 5.0, and 5.1. Fig- 
ure [4] clarifies that for b < 5crit ~ (5.0 ± 0.1)moT/v, a 
BH is formed after the collision for the chosen velocity, 
0.6 < V < 0.9. Extrapolating the result for w ^ 1 under 
the assumption that the discovered relation holds even for 
u — !■ 1 , the maximum impact parameter is determined to 
be (2.50±0.05)A/o. 

Figure plots outgoing part of the Newman-Penrose 
quantity for v = 0.8 and bv/rrioT = 5.088. I = m = 2 and 
I = m = A modes are plotted for the best-resolved run 
with h = 0.042mo. For I — m = 2, we also show the re- 
sult for h — 0.05r7io, which agrees with the best-resolved 
result with a small error. Note that the waveforms are 
qualitatively similar irrespective of the value of v as far 
as 6 w fcciit, although the maximum amplitude steeply 
increases when b approaches bcrit ■ Figure shows that 
gravitational waves are efficiently emitted after the onset 
of collision: When two BHs approach each other, am- 
plitude of gravitational waves gradually increases. As 
the separation of two BHs becomes sufficiently small, 
they constitute a bound orbit and quasiperiodic gravi- 
tational waves are emitted. After substantial fraction of 
gravitational waves are emitted, two BHs merge to be 
a new BH, and then, ring-down gravitational waves as- 
sociated with fundamental quasinormal modes are emit- 
ted for tj-ct ^ 40Afo. In this case, the formed BH is 
rapidly rotating with the spin parameter « 0.73 ± 0.02, 
and hence, the damping time scale is longer than that 
for nonrotating BHs [12]. Because the orbital velocity is 
very large, higher-multipole components of gravitational 
waves are also enhanced significantly (cf. the waveform 
for I = m = 4:). 

Total energy and angular momentum dissipated by 
gravitational waves are « 25 ± 5% and « 65 ± 5% of 
the initial energy and angular momentum, respectively, 
for b ^ 6crit and for v — 0.9. The totally emitted grav- 
itational radiation for b ~ 5crit slightly decreases with 
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decreasing but still, AE/Mq > 20% and A J/ J > 60% 
for 0.6 < V < 0.8. I — \m\ — 4 modes contribute to 
AE and AJ by w 10-15% and by w 15-20%, respec- 
tively, for b ^ 6crit- It should be noted that in the limit 
b —f &crit, total amount of gravitational radiation may 
be slightly larger than that presented here, because the 
lifetime of the formed binary orbit could be longer. We 
note that the results for the total amount of gravita- 
tional radiation are consistent with the mass and spin of 
the BH finally formed within an acceptable error for the 
best-resolved runs: The mass and angular momentum 
of the BHs estimated from apparent horizon are always 
smaller than those expected from gravitational radiation 
(i.e., Mo - AE and J - AJ) by < 0.05Mo and < 0.1 J, 
respectively, for the best-resolved run. The error is larger 
for the larger value of v. The reason for this error is that 
energy and angular momentum are dissipated spuriously 
by numerical effects associated with finite grid resolu- 
tion. However, our results show a behavior of conver- 
gence slightly better than second order with improving 
grid resolution. 

IV Discussion and summary: We find that the 
largest value of the impact parameter for the BH for- 
mation after the collision is 6crit ~ (2.50 ± 0.05)Mo/w. 
For such value of the impact parameter, the initial value 
of the spin parameter of the system is 



J 

Ml 



1.25 ±0.03. 



(18) 



For the BH formation, the spin parameter of the formed 
BH should be smaller than unity if the cosmic censor- 
ship holds [l3|. This implies that a large fraction of 



the angular momentum is dissipated by gravitational 
radiation during the collision. We estimate the total 
amount of the angular momentum and energy dissipated 
by gravitational radiation are AJ — (0.65 ± 0.05) J and 
AE — (0.25 ± 0.05)Afo, respectively. The expected spin 
parameter of the formed BH is 



J - AJ 
{Mo - AE)^ 



(0.6 ±0.1)^. 



(19) 



Thus, even if the spin parameter of the system is initially 
1.25, the resulting value at BH formation is smaller than 
unity. As this discussion clarifies, gravitational radiation 
increases the impact factor by ~ 25%, and as a result, 
the cross section by ^ 50%. It is also worth noting that 
the formed BH does not have an extremely large spin 
^ 1, but approximately 0.8 ± 0.1 even for b < 6crit- 

In this work, we adopt initial data which satisfies the 
constraint equations only approximately. Although the 
violation is tiny (cf. Fig. [1]), this error produces a small 
error in estimation of the critical cross section, and AE 
and A J of gravitational waves. To determine these quan- 
tities strictly, it is necessary to perform simulations using 
improved initial condition. 

This work is a first step toward detailed understand- 
ing of high-velocity collision of two BHs in higher- 
dimensional spacetime. We plan to develop a numerical 
code for higher-dimensional spacetimes. 
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